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Abstract. In this paper we propose a new algorithm for obtaining the rational integrals of 
the full Kostant-Toda lattice. This new approach is based on a reduction of a bi-Hamiltonian 
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The Toda lattice is arguably the most fundamental and basic of all finite dimensional integrablc 
systems. It has various intriguing connections with other parts of mathematics and physics. 
The Hamiltonian of the Toda lattice is given by 



N N-l 

H( qi , ...,q N , Pl ,... ,p N ) = ]T -p? + ^ e9 '~ qi+1 - (1) 



=1 



in 
o 

•l-H 
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| Equation Q is known as the classical, finite, non-periodic Toda lattice to distinguish the 

^ ' system from the many and various other versions, e.g., the relativistic, quantum, periodic etc. 

This system was investigated in [^1 H3 El OH EH 0^] and numerous of other papers; see [1] 
for a more extensive bibliography. 
Hamilton's equations become 

pj = e <U-i-<U _ e ©-©+i. 

The system is integrable. One can find a set of independent functions {-Hi, . . . , H^} which are 
constants of motion for Hamilton's equations. To determine the constants of motion, one uses 
Flaschka's transformation: 

o i = l e 3(*-*+i) ) b i = - V -p l . (2) 

Then 

di = ai(b i+ i - bi), bi = 2(aj - a|_i). (3) 
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These equations can be written as a Lax pair L = [B,L], where L is the Jacobi matrix 
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(4) 



and B is the skew-symmetric part of L. This is an example of an isospectral deformation; the 
entries of L vary over time but the eigenvalues remain constant. It follows that the functions 
Hi = j tr L l are constants of motion. 

Note that the Lax pair (jlj) has the form 



L(t) = [PL(t),L(t)}, 



where P denotes the projection onto the skew-symmetric part in the decomposition of L into 
skew-symmetric plus lower triangular. 

The Toda lattice was generalized in several directions. 

We mention the Bogoyavlensky-Toda lattices which generalize the Toda lattice (which cor- 
responds to a root system of type A n ) to other simple Lie groups. This generalization is due 
to Bogoyavlensky pp. These systems were studied extensively in ^S] where the solution of the 
systems was connected intimately with the representation theory of simple Lie groups. 

Another generalization is due to Deift, Li, Nanda and Tomei |5j who showed that the system 
remains integrable when L is replaced by a full (generic) symmetric n x n matrix. 

Another variation is the full Kostant-Toda lattice (FKT) 121]. We briefly describe the 

system: In Kostant conjugates the matrix L in (jlj) by a diagonal matrix to obtain a matrix 
of the form 



(h 


1 





ai 


b-2 


1 





a 2 


h 



: '•■ '•■ ■•• 

I '•■ '•• 1 

\0 a n _i b n ) 

The equations take the form 
X(t) = [X(t),PX(t)], 



where P is the projection onto the strictly lower triangular part of X(t). This form is convenient 
in applying Lie theoretic techniques to describe the system. 

To obtain the full Kostant-Toda lattice we fill the lower triangular part of L in (jlj) with 
additional variables. (P is again the projection onto the strictly lower part of X(t)). So, using 
the notation from an d HH 



X(t) = [X(t),PX(t)}, 
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where X is in e + i?_ and PX is in iV_. _B_ is the Lie algebra of lower triangular matrices 
and iV_ is the Lie algebra of strictly lower triangular matrices. The fixed matrix e has a general 
form in terms of root systems: 



^ ] X Oti 



where A denotes the set of simple roots. In fact the FKT lattice itself can easily be generalized 
for each simple Lie group; see [2J Ell- 
in the case of sl(4, C) the matrix X has the form 



X 



(fx 


1 





°\ 


91 


h 


1 





hi 


92 


h 
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\h 


h 2 


93 


fj 



(6) 



with y: fi = o. 

i 

The functions Hi = I Tr X i are still in involution but they are not enough to ensure integra- 
bility. This is a crucial point: the existence of a Lax pair does not guarantee integrability. There 
are, however, additional integrals which are rational functions of the entries of X. The method 
used to obtain these additional integrals is called chopping and was used originally in for the 
full symmetric Toda and later in [H] for the case of the full Kostant-Toda lattice. 

In this paper we use a different method of obtaining these rational integrals which does 
not involve chopping. This method uses a reduction of a bi-Hamiltonian system on gl(N, R), 
a system which was first defined in |2U1 121j . In A. Meucci presents the bi-Hamiltonian 
structure of Toda3, a dynamical system studied by Kupershmidt in (TJj as a reduction of the 
KP hierarchy. Meucci derives this structure by a suitable restriction of the set of maps from Zj, 
where Zj is the cyclic group of order d, to GL(3,R), in the context of Lie algebroids. 

In (20] the bi-Hamiltonian structure of the periodic Toda lattice is investigated by the reduc- 
tion process described above using maps from Z& to GL(2,M). This approach parallels the work 
of [Zj where the continuous analog of the Toda lattice is studied, namely the KdV. If instead the 
target space is gl(3,M) one obtains the Boussineq hierarchy. The work of Meucci is a discrete 
version of this approach. If one generalizes the cases N = 2, 3, i.e. consider maps from Z& to 
GL(N,M) the resulting system will be denoted by TodaAr. In the present paper we use the 
results of |2()[ 121) as a starting point. We begin with the bi-Hamiltonian system obtained on 
gl(N, K) in the particular case d = N and use a further reduction to obtain the dynamics of the 
full Kostant-Toda lattice. We then propose a new algorithm which produces all the rational 
integrals for the FKT lattice without using chopping. We present this algorithm by specific 
examples (N = 3, 4) in Section 4. Sections 2 and 3 contain a general review of the old methods 
and results on integrability and bi-Hamiltonian structure of the FKT lattice. 



2 Integrability of the FKT lattice 

Let Q = sl(n), the Lie algebra of n x n matrices of trace zero. Using the decomposition 
Q = B + we can identify Bl. with the annihilator of iV_ with respect to the Killing form. 
This annihilator is Thus we can identify Bl. with i?_ and therefore with e + £?_ as well. 
The Lie-Poisson bracket in the case of sZ(4, C) is given by the following defining relations: 



{9i,9i+i} = hi, 
{9i,fi} = ~9i, 
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Wi, 


t \ 

H+lf 


= 9i, 


{hi, 


hi = 


-hi, 


{hi, 


t \ 

Ji+2} 


= hi, 


r 

\9i, 


M = 


fa, 


{93, 


hi} = 


-fa, 


{fa 


fi} = 


-fa, 


{fa 


h} = 


fa. 



All other brackets are zero. Actually, we calculated the brackets on gZ(4, C); the trace of X now 
becomes a Casimir. The Hamiltonian in this bracket is H2 = ^ Tr X 2 . 

Remark 1. If we use a more conventional notation for the matrix X, i.e. Xij for i > j, xa+i = 1, 
and all other entries zero, then the bracket is simply 

{xij, x k i} = 5iiX k j - 5j k Xii. (7) 



The functions Hi = i Tr X % are still in involution but they are not enough to ensure integra- 
bility. There are, however, additional integrals and the interesting feature of this system is that 
the additional integrals turn out to be rational functions of the entries of X. We describe the 
constants of motion following references [f3 I24j . 

denote by {X — A Id) a.) the result of removing the first k rows and 



For k = 0, . . . 



(n-l) 
2 



last k columns from X — A Id, and let 

det(X - AId) (fc) = E 0k \ n - 2k + ■■■ + E n _ 2kM . 
Set 

det(X-AId) (fc) n _ 2fc x„_2fc-i, , j 

^ — A + \~ 1 n-2k,k- 

The functions I rk , r = 1, . . . , n — 2k, are constants of motion for the FKT lattice. 
Example 1. We consider in detail the gl(3, C) case: 

X = 

Taking H2 = \ tr X 2 as the Hamiltonian, and the above Poisson bracket 

x = {H 2 ,x} 
gives the following equations: 

A = -9i, 
h = 91 ~92, 
h = 92, 

91 = 9l{h ~ h) ~ h i, 

92 = 92U2 ~ f-i) + h, 
h^hiih-fs) 
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Note that H\ = f\ + f 2 + fz while H 2 = 5 (/1 + /| + /f) + Si + 92- These equations can be 
written in Lax pair form, X = [H,X], by taking 




B 



The chopped matrix is given by 

9i /2 - A 
hi 92 

The determinant of this matrix is h\X + gig 2 — ^1/2 and we obtain the rational integral 

T 9i92 ~ hif 2 , Q s 
111 = h{ • (8) 

Note that the phase space is six dimensional, we have two Casimirs (Hi, In) and the functions 
(H 2 ,Hz) are enough to ensure integrability. 

Example 2. In the case of gl(4,C) the additional integral is 

T 919293 ~ 31/3/12 - f293h + hih 2 , , , 

hi = ; h /2J3 - 92- 

fci 

and 

T gih 2 + g3hi . , 

1 n — ; J2 — J3 

is a Casimir. 

In this example the phase space is ten dimensional, we have two Casimirs (Hi, In) and the 
functions (H 2 , H3, H4, /21) are independent and pairwise in involution. 



3 Bi-Hamiltonian structure 

We recall the definition and basic properties of master symmetries following Fuchssteiner |13j . 
Consider a differential equation on a manifold M defined by a vector field \- We are mostly 
interested in the case where x is a Hamiltonian vector field. A vector field Z is a symmetry of 
the equation if 

[z,x]=o. 

A vector field Z is called a master symmetry if 

[[z,x\,x\=o, 

but 

Master symmetries were first introduced by Fokas and Fuchssteiner in |12j in connection with 
the Benjamin-Ono Equation. 

A bi-Hamiltonian system is defined by specifying two Hamiltonian functions Hi, H 2 and 
two Poisson tensors tt\ and tt 2 , that give rise to the same Hamiltonian equations. Namely, 
iti'VH 2 = n 2 VH±. The notion of bi-Hamiltonian system was introduced in ^S] in 1978. 
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Another idea that will be useful is the Gelfand-Zakharevich scheme for a pencil of Poisson 
tensors. Consider a bi-Hamiltonian system given by two compatible Poisson tensors 7Q, 7r 2 . 
Compatible means that the pencil 

{ , }a = { , }tu + A{ , } 7T2 

is Poisson for each value of the real parameter A. Gelfand and Zakharevich in Jl] consider the 
special case in which the pencil possesses only one Casimir. Under some mild conditions they 
prove the following: 
Let 

F x = F + XF l + --- + X n F n . 

Then Fq is a Casimir for tt\, F n is a Casimir for 7r 2 and, in addition, the functions Fq, F\, . . . , F n 
are in involution with respect to both brackets tti and 7T2- 

We now return to the FKT lattice. We want to define a second bracket 1x2 so that H\ is the 
Hamiltonian and 

7T 2 V-ffl = 7T1 VH 2 . 

i.e. we want to construct a bi-Hamiltonian pair. We will achieve this by finding a master 
symmetry X\ so that 

Xi(TrA^) =iTvX i+1 . 
We construct X\ by considering the equation 

X = [Y,X]+X 2 . (9) 
We choose Y in such a way that the equation is consistent. One solution is 

n n—1 



Y = (XiEii + & E i,i+l> 



i=l i=l 

i-1 

where A =i, ai = ifo + Yl fk- 

k=l 

The vector field X\ is defined by the right hand side of @. 
For example, in 5/(4, C) the components of X\ are: 

X l (f l ) = 2g 1 + fl 

X 1 (f 2 ) = 3g 2 + fl 

Xi(h) = -g2 + 4g 3 + fl, 

X 1 (f i ) = -2g 3 + fi, 

Xi(gi) = 3hi + gifx + 351/2, 

Xi(g 2 ) =4/i2 + 4 52 /3, 

-^1(33) = -h2 - 53/3 + 553/4, 

Xi(hi) = 5152 + 4fci + hifi + hif 2 + 4/11/3, 

Xi(h 2 ) = 5253 + ^2/3 + 5/12/4, 

^l(fcl) = 53^1 + 51 ^2 + fcl/l + hf2 + fcl/3 + 5&i/4. 
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The second bracket 112 is denned by taking the Lie derivative of tt\ in the direction of X\. 
The bracket 712 is at most quadratic, i.e. in the case n = 4 



r ^ 

\9i,9i+it 


1 U -P 

— 9i9i+l + 


{9i,hi} = 


7 

gm, 


{g i+1 ,hi} 


7 

— -9i+l n i, 


{9i,fi} = 


~9ilii 


{gufi+i} 


= 9ifi+l, 


{9i,fi+2} 


= hi, 


{9i+i,fi} 


= -hi, 


{hi,fi} = 


hi fii 


{hi, fi+2} 


= hifi + 2, 


{fii fi+l} 


= 9i- 



This bracket was first obtained in [3] using the method described above (i.e. master symme- 
tries). A closed expression for this bracket was obtained later by Faybusovich and Gekhman 
in [S]. It was obtained using i?-matrices and the expression takes the following simple form: 



sign(/c - i)xijx k j if j = I, 

%ij%il if ^ — ^; 

x%jXki -\- xnXkj if i < k < j., 

x u if k = j + 1. 



(10) 



As we will see in the next section, there is a linear and a quadratic bracket on gl(n,M) whose 
restriction to the FKT lattice coincides with the brackets tt\, tt2- 



4 A new approach 

We already mentioned in the introduction that our starting point is the work of A. Meucci |2()1 
I21j . His work is a discrete analogue of a procedure that produces the KdV, Boussineq and 
Gelfand-Dickey hierarcies 0. For example, the KdV is bi-Hamiltonian and it can be obtained 
by reducing the space of C°° maps from S 1 to gl(2, M). If one considers the space of maps 
from S 1 to gl(3, M) the Boussineq hierarchy is obtained. In |2()l 121 j the discrete version of the 
procedure is considered. The circle, S , is replaced by the cyclic group and one considers 
maps from to GL(N, R). One obtains, after reduction, equations that are bi-Hamiltonian. In 
the case N = 2 the resulting system is the periodic Toda lattice and for iV = 3 a system studied 
by Kupershmidt in J7j . We will not get into the details of the procedure but rather we will use 
the results as our starting point for our own purposes. We will content with a short outline of 
the constructions of [20] and The basic object is the space of maps from to GL(N,M). 
This space is endowed with a structure consisting of a Poisson manifold together with a pair 
of Lie-algebroids suitably related. The next step is a Marsden-Ratiu type of reduction and the 
result is a bi-Hamiltonian system with a pair of Poisson structures on gl(N,M). The Lax pair 
of the resulting system contains two spectral parameters and the theory of Gelfand-Zakarevich 
applies. The system turns out to be integrable with the required number of integrals. We give 
explicit formulas that we have computed from j2Hi m the case N = 3 both for the pair of Poisson 
brackets, the Lax pair and the integrals of motion. In the case N = 4 we display the Lax pair 
and the polynomial integrals of motion. To obtain the FKT lattice one has to perform a further 
reduction to the phase space of the FKT lattice and to obtain the rational integrals we propose 
a new algorithm which is the central new result of our paper. At the present we do not have 
a Lie algebraic interpretation of this algorithm. We illustrate with two examples: 



8 



P.A. Damianou and F. Magri 



Example 3 (The gl(3,C) case). The phase space of the system obtained by the procedure 
of Meucci is nine dimensional, i.e. matrices of the form 



We compute the pair of Poisson brackets on the extended space with variables 

{fl, f2, /3j 9i, 92, 93, h±,h 2 , h 3 }. 
The Lie-Poisson bracket is defined by the following structure matrix 














-9i 





.93 


-hi 


h 2 


^ 













.91 


~92 








-h 2 


h 3 
















92 


-93 


hi 





-h 3 




.9i 


-9i 








-hi 


h 3 
















.92 


~92 


hi 





-h 2 













-93 





93 


-h 3 


h 2 
















hi 





-hi 






















-h 2 


ho 




























-h 3 


h3 

















o J 



and the quadratic Poisson bracket is defined by A — A 1 where A is the matrix 





91 


-93 


~9i fl 


hi - h 2 


93 fl 


-hih 


h 2 fi 


> 








92 


91 f 2 


-92h 


h 2 - h 3 





-h 2 f 2 


/13/2 











h 3 - hi 


92 J 3 


-53/3 


hih 





-h 3 f 3 














-hif 2 - gig 2 


h 3 fi +gig3 


-gihi 





9ih 3 

















-h 2 f 3 - g 2 g 3 


92h 


-.92/12 


























93h 2 


-53/13 
























































\o 























J 



The Lax matrix with two spectral parameters is given by 





Let p(X, fi) 



detL^x. Write 

-fJ 9 + c 2 (X)fi 6 + ci(X)fi 3 + c (X). 



p{X,fi) 



Then 



02(A) = A 3 + k 2 X 2 + kiX + fco, 



where 



k 2 
ki 
ko 



fi + h + h, 

/1/2 + fih + f2h + 91+92 + 93, 

/1/2/3 + fi92 + f293 + h9l +hi + h 2 + h 3 



Full Kostant-Toda Lattice 



9 



and 

ci(A) = h\ + / , 

where 

h = ~h 3 g 2 - g 3 hi - h 2 gi, 

k = -fi92h 3 + 515352 - 51^2/3 - hg 3 f 2 - h 2 h 3 - h 3 hi - h x h 2 . 
Finally, 

c (A) = h^hs. 

The functions ki, li, Co are all in involution in the Lie-Poisson bracket. The functions k 2 , h, 
Co are all Casimirs. 

In the quadratic bracket ki, li, and Co are all in involution, ko, Iq and Co are Casimirs. 
Let 

l _ k _ -fi92h 3 + giff352 ~ 51^2/3 ~ ^153/2 ~ h 2 h 3 ~ fahi - h x h 2 
h -h 3 g 2 - g 3 hi - h 2 g x 

Setting h 2 = h 3 = in I we obtain In (JSJ). 

Example 4 (The gl(4,C) case). In order to give a complete comparison of the previous 
and the present method of obtaining the rational invariants we consider first integrability using 
chopping. 

We consider the matrix L given by 



L 



(h -10 0\ 

91 h -1 

hi 92 f 3 -1 

\h h 2 g 3 Uj 



Note that we are using — e instead of e in order to get the integrals to match exactly. 
In this case the chopped matrix has the form 

1 9i /a -A -1 
Chi(X)=\hi 52 /3-A 
\h h 2 g 3 

The characteristic polynomial has the form 

hX 2 + (gih 2 + hig 3 - kif 2 - hf 3 )X + 515352 - 51^2/3 - hif 2 g 3 - h x h 2 + fci/2/3 + hg 2 . 
We obtain the following two rational invariants 
h 2 gi + 9zhi 



and 



"251 -1- 53 (f , f 
111 = 7T (72 + h 



5i5352 - 5i^2/3-/ii/253 -/ii/i2 . , . . 
*2l = ; V 7273 +52- 
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We now turn to the gentle approach, i.e. integrability without chopping. Consider the fol- 
lowing Lax pair with two spectral parameters: 

/ (/i + A)/i 3 k 2 - // ksn 

T giV 2 (/2 + A)/x 3 h-fi 4 

h lf x g 2 f (/ 3 + A)/i 3 

\ h\ - // h 2 fj, g 3 n 2 

Taking determinant we obtain the polynomial 

P\,^ = V 6 + K 3 {X)^ 12 + K 2 (X)fi 8 + Ki(A)/i 4 + Kq(X). 
We present the explicit expressions for the polynomials -?Q(A). 

. K 3 (X) = K 33 X 3 + K 32 X 2 + K 31 X + K 30 , 

where 

^33 =h+h+h+ u, 

K 32 = gi + 92 + 93 + 94 + /1/2 + hh + /1/4 + f-ih + hf± + hh, 

K 31 =h 1+ h 2 + h 3 + h 4 + hhh + /l/2/4 + /l/3/4 + /2/3/4 + /294 

+ /394 + hgi + hgi + /452 + /253 + fm + fm, 

K 30 = h + k 2 + k 3 + k 4 + /1/2/3/4 + 5if3 + 5254 + M3 + Ml 

+ M2 + /3/14 + /2/354 + /3/451 + /l/452 + /l/253- 

• K 2 (X) = K 22 X 2 + K 21 X + K 20 , 

where 

#22 = -h 2 h± - 52^3 - 53^4 - hh 3 - gik 2 - hg A , 

K 21 = g 2 h A g 3 - gik 2 f 4 - fih 2 h A - kig A f 3 - fci/254 - hg 2 k 3 - h 2 h A f 3 - g±k 2 f 3 
+ 5152/13 + hgm + 51^254 - Zi53&4 - hf 2 h 3 - g 2 k 3 fi - f 2 g 3 k A - h 2 k 4 

- h 2 k 3 — hih 3 f4 — k 2 h>4 — h 3 k A — h\k 2 — h\k 3 — k\h 3 — k\h A , 
K20 = -k 2 k A - k 2 k 3 - k±k 2 - k\k 3 - kik 4 - k 3 k 4 - kif 2 g A f 3 - f\h 2 h A f 3 

+ /152/1453 - hik 2 U ~ fmhfi - fif 2 g 3 k A - gig 2 g 4 g 3 + g\h 2 g 4 f 3 

- hif 2 h 3 f A + g\g 2 h 3 f A + hif 2 g A g 3 - g\k 2 f 3 f A - f\h 2 k 3 - g\k 2 g 3 - f\h 2 k A 

- k\g 2 g A + hih 2 g A - k 2 h A f 3 - f 2 h 3 k A + g 2 h 3 h A - g 2 g A k 3 - g\g 3 k A + g±h 2 h 3 

- hik 3 f A + hih A g 3 - kif 2 h 3 - hh A f 3 . 

. K 1 (X) = K 11 X + K W , 

where 

Ku = h\k 2 k 3 + h 2 k 3 k A + k\k 2 h A + k\h 3 k A , 

K 10 = hik 2 k 3 f A + fih 2 k 3 k A - k\g 2 h 3 h A + hg 2 g A k 3 + k x k 2 h A f 3 + fci/2/13^4 

- hxh 2 g A k 3 + hih 2 h 3 h A - hik 2 h A g 3 - g\h 2 h 3 k A + g\k 2 g 3 k A + k x k 2 k A 
+ k\k 2 k 3 + k 2 k 3 k 4 + kik 3 k A . 

• K (X) = -k\k 2 k 3 k A . 



54/T \ 
h A fj, 
k A — 
(/4 + AV 3 / 



(11) 
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Remark 2. We note that K 3 (X) gives the polynomial invariants. Clearly H\ = K 33 = tr L. We 
also have H 2 = \txL 2 = ±Kf 3 - K 32 and H 3 = ± tr L 3 = ±iff 3 - K 33 K 32 + K 31 . 
Finally H 4 = + K 33 K 31 + \K 32 - K 30 - K 33 K 32 . 

These last relations hold provided that k 2 = k 3 = k^ = 0, h 3 = /14 = and 54 = 0. 

Remark 3. The next coefficient K 2 (\) gives the rational invariants. We form the quotient 
and set k 2 = k 3 = k& = and h 3 = /14 = 0. We obtain 

K 2 \ -54(^1/3 + ^1/2 - h x g 3 - gih 2 ) {h x g 3 + g±h 2 ) 

~zF~ ~ ; - J2 + h ; ■ 

A 22 -kigi h 

This is precisely —in. 

Similarly we form to obtain precisely i 2 i- 

Remark 4. The last two terms, namely K\(X) and Kq(X) become identically zero once we set 
k 2 = k 3 = /c4 = 0, h 3 = /14 = 0. 

In the general case, the polynomial involves polynomials j4(A), i?i(A), . . . , Bfr(X) and 

Ci(A), . . . , C s {\) with k = 1 ^ and s = n — k — 1. The polynomial A(X) can be used to 

obtain the polynomial integrals. The polynomials B{(\) give the rational integrals using the 
procedure described above and the Cj(X) vanish identically once we restrict to the phase space 
of the FKT lattice. A detailed proof of the general case will be given in a future publication. 
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